Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


[J.  Res.  Natl.  Inst.  Stand.  Technol.  99,  475  (1994)] 


A Random  Field  Excursion  Model  of 
Salt  Induced  Concrete  Delamination 


Volume  99 


Number  4 


July- August  1994 


Gordon  A.  Fenton 

Technical  University  of  Nova 
Scotia, 

Halifax,  Nova  Scotia,  B3J  2X4 


Salt  induced  concrete  delamination  is  a 
problem  often  encountered  in  parking 
garage  slabs  in  northern  climates  where 
deicing  salts  are  heavily  used.  A ran- 
dom field  model  of  the  delamination 
process  is  investigated  and  compared 
against  some  field  data.  Delaminated 
regions  of  the  slab  are  modeled  as  ex- 
cursions of  a random  field  above  a pre- 
scribed threshold  and  the  growth  of 
these  regions  with  time  is  obtained  by 
allowing  the  threshold  to  fall  as  a func- 


tion of  time.  Simulation  based  excur- 
sion statistics  are  used  to  obtain  the 
mean  and  variability  of  various  aspects 
of  the  delamination  process  using  this 
model. 
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1.  Introduction 

The  delamination  and  spalling  of  concrete  sur- 
faces in  parking  structures  in  northern  climates  is 
an  ongoing  and  expensive  problem.  The  Trans- 
portation Research  Board  [1]  estimates  that  be- 
tween 50  and  150  parking  structures  in  the 
Northeast  and  Midwest  United  States  will  need  to 
be  rehabilitated  each  year  for  the  next  10  years  at 
an  average  cost  of  $1  million  per  structure.  With- 
out proper  design  and/or  maintenance,  deicing 
salts  brought  in  by  vehicles  from  the  roadways  are 
deposited  on  the  concrete  surface  along  with  water. 
Chloride  ions  gradually  penetrate  the  concrete  and 
electro-chemical  processes  lead  to  corrosion  of  the 
reinforcing  steel.  This  results  in  both  a degradation 
of  structural  integrity  and,  since  the  corrosion 
products  occupy  considerably  more  volume  than 
the  original  steel,  delamination  or  spalling  of  the 
concrete  surface  and  loss  of  utility. 

This  paper  is  aimed  at  developing  a tool  to  aid 
in  a rational  probabilistic  approach  to  the  rehabil- 
itation of  parking  structures.  Such  an  approach 
allows  the  optimal  allocation  of  limited  resources 


to  this  ongoing  and  rather  expensive  maintenance 
issue.  A simple  stochastic  model  involving  only  a 
few  parameters  is  used  herein  to  represent  what  is 
known  to  be  a complex  phenomenon.  The  follow- 
ing factors  are  suggested  by  Public  Works  Canada 
[2]  to  have  the  largest  effect  on  the  onset  of  con- 
crete delamination; 

• chloride  ion  input:  quantities  of  deicing  salts 
used, 

• concrete  permeability:  influenced  in  turn  by 
water/cement  ratio,  intensity  and  frequency  of 
cracking,  surface  coatings/sealers  and  construc- 
tion practices, 

• ambient  temperature,  humidity,  precipitation 

• concrete  cover  depth 

• pH  of  local  aggregates 

• conductivity:  wetness  of  concrete. 

It  is  immediately  recognized  that  virtually  all  of 
these  factors  are  highly  variable  from  structure  to 
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structure  and  within  a single  structure  from  point  to 
point.  As  well,  on  a practical  basis  for  an  existing 
structure,  some  of  the  factors  are  unknowable  ex- 
cept through  extensive  destructive  testing.  Clearly  a 
representative  model  should  not  require  experi- 
mental validation  on  a per  structure  basis,  where 
the  experimentation  may  be  more  expensive  than 
the  final  repairs.  This,  in  fact,  is  the  primary  motiva- 
tion for  the  use  of  stochastic  models. 

In  a different  approach  to  the  same  problem, 
Attwood  et  al.  [3]  develop  a limit  state  function  for 
parking  structures  based  on  a critical  fraction  of 
floor  area  delamination,  Dcr, 

m =Dcr- aS{t -t\),  t>ti,  (1) 

where  m <0  denotes  the  “failure”  state,  S is  the 
annual  delamination  rate,  a is  a correction  factor 
accounting  for  ambient  temperature  and  crack 
widths,  t is  time  in  years,  and  h is  the  time  to  initia- 
tion of  delamination  (also  in  years).  Note  that  the 
second  term  is  taken  to  be  zero  for  all  t<ti. 
Attwood  et  al.  employ  a First-Order  Reliability 
Method  (FORM)  to  estimate  reliabilities  associ- 
ated with  the  delamination  of  a parking  garage 
structure.  In  their  approach,  all  of  the  factors  ap- 
pearing in  Eq.  (1)  are  expressed  in  terms  of  random 
variables  having  assumed  distributions  and  the  joint 
cumulative  probabilities  are  evaluated  using 
FORM. 

In  this  paper  a 2-dimensional  random  field  model 
representing  the  spatial  delamination  process  over 
time  is  investigated.  The  random  field  may  be 
loosely  interpreted  as  the  out-of-plane  stress  field 
at  the  reinforcement  level  (which  changes  randomly 
from  point  to  point  over  the  area  of  the  slab).  De- 
laminated regions  are  represented  by  excursions  of 
the  random  field  above  some  threshold  which  can 
be  thought  of  as  a critical  concrete  resistance  to 
horizontal  splitting.  In  fact  the  stress  field  itself  will 
not  be  considered  in  that  it  is  unmeasurable  in  prac- 
tice. Rather,  the  field  excursions  will  be  used  to  rep- 
resent the  delaminated  regions  directly,  since  these 
regions  are  measurable  to  some  extent.  A simple 
random  field  model  is  adopted  whose  primary  moti- 
vation is  to  attempt  to  shed  light  on  the  following 
questions; 

1)  what  is  the  mean  and  variability  of  the  total 
area  of  delamination  as  a function  of  time? 

2)  what  is  the  average  size  of  individual  delami- 
nated regions  as  a function  of  time? 

3)  how  many  delaminated  regions  can  be 
expected  in  a slab? 


These  are  essentially  questions  regarding  the  statis- 
tics of  excursion  regions  and  so  some  simulation 
results  regarding  excursion  statistics  will  be  pre- 
sented in  the  next  section. 

Figure  1 illustrates  an  example  of  excursions  of  a 
random  field  above  some  predefined  threshold.  In 
the  context  of  this  paper,  the  dark  regions  can  be 
viewed  as  areas  where  concrete  delamination  has 
occurred  at  some  fixed  time.  As  time  progresses, 
delaminated  regions  are  expected  to  grow  in  size, 
corresponding  to  a falling  threshold  level.  A falling 
threshold  is  equivalent  to  a rising  mean,  in  the  case 
of  a homogeneous  field.  Since  the  excursion  statis- 
tics in  the  next  section  are  developed  as  functions  of 
the  threshold,  the  falling  threshold  interpretation  is 
used  here  rather  than  a rising  mean.  In  either  case, 
a nonhomogeneous  field  could  be  employed,  if  the 
data  so  indicated,  by  considerably  extending  the 
simulation  based  study  (in  the  absence  of  analytical 
results).  In  this  preliminary  investigation,  only  a ho- 
mogeneous field  is  considered. 


Fig.  1.  Excursions  of  a two-dimensional  random  field  above 
threshold  ba. 

To  represent  the  delamination  process,  excur- 
sions of  an  isotropic  Gaussian  random  field  having 
mean  zero  and  unit  variance  will  be  used.  Although 
other  values  of  mean  and  variance  are  possible,  the 
excursion  statistics  are  dependent  purely  on  the 
distance  between  the  mean  and  the  threshold.  This 
distance  can  be  expressed  in  units  of  a,  the  stan- 
dard deviation,  so  that  there  is  no  advantage  in 
choosing  anything  other  than  a mean  zero,  unit 
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variance  field.  The  quantity  ba  will  be  referred  to 
henceforth  as  the  physical  threshold  and  b alone  as 
the  threshold.  The  choice  of  an  isotropic  Gaussian 
process  has  been  made  largely  for  simplicity,  there 
being  little  evidence  available  to  clearly  justify  other 
types  of  random  functions. 

2.  Excursion  Statistics  in  Two  Dimen- 
sions 


between  cell  variance  and  point  variance  will  be  ig- 
nored in  this  paper,  although  the  results  presented 
in  the  various  plots  to  follow  are  accurate  in  this  re- 
spect. For  a more  rigorous  treatment  of  this  issue, 
see  Ref.  [6]. 

Within  a given  domain  V — [0,  Li]  x [0,  L2]  of  area 
At=Li  L2,  the  total  excursion  area  per  unit  area, 
Db,  where  the  process  Z(jc)  exceeds  some 
threshold,  can  be  defined  by 


With  respect  to  excursion  statistics,  such  as  the 
mean  number  and  area  of  isolated  excursions,  ana- 
lytical results  developed  to  date  are  asymptotic  in 
nature,  accurate  only  at  very  high  thresholds  where 
the  excursion  process  approaches  a Poisson  point 
process.  Often  in  engineering  problems  the  interest 
is  in  thresholds  which  are  quite  a bit  lower,  such  as 
the  delamination  process  considered  herein.  In  or- 
der for  the  proposed  random  field  model  to  be  use- 
ful in  this  context,  excursion  statistics  should  be 
available.  This  section  summarizes  a study  in  which 
excursion  statistics  are  obtained  through  Monte 
Carlo  simulation.  Specifically,  an  ensemble  of  2000 
realizations  of  an  isotropic  zero  mean,  unit  variance 
Gaussian  random  field,  Z(jc),  with  Markovian  co- 
variance  function, 

R(T)  = o-V(T)  = o-^expj-||T|j,  (2) 


Db=^  f I {Z(x)-ba)  dx, 
Jy 


(3) 


where  ba  is  the  physical  threshold  of  interest,  a^  be- 
ing the  variance  of  the  process,  and  /(•)  is  the  indi- 
cator function  defined  on  V (taken  to  be  zero 
outside  the  domain  V) 


if  r>0 
if  r < 0 • 


(4) 


For  a homogeneous  process,  the  expected  value  of 
Db  is  simply 

E[Db'\  = P[Z>bal  (5) 

which,  for  a zero-mean  Gaussian  process  yields 


E[Db]=\-nb),  (6) 


are  produced  using  the  Local  Average  Subdivision 
(LAS)  method  [4,  5],  where  t is  the  lag  vector  and 
d is  the  scale  of  fluctuation.  The  scale  of  fluctuation 
is  loosely  interpreted  as  the  distance  over  which 
correlation  is  significant.  Since  many  of  the  statis- 
tics of  interest  depend  strongly  on  the  scale  of  fluc- 
tuation, 2000  realizations  were  generated  at  each  of 
5 different  scales  of  fluctuation. 

Individual  realizations  are  decomposed  into 
excursion  regions  and  “holes,”  using  a space-filling 
algorithm,  over  a range  of  thresholds  b = [-4,  4]. 
The  mean  and  variance  statistics  of  the  excursions 
are  estimated  over  the  ensemble.  It  should  again  be 
emphasized  that  b is  measured  in  units  of  standard 
deviation,  for  example  b = 2 implies  a threshold  at 
two  standard  deviations  from  the  mean.  For  a unit 
variance  field  there  is  no  distinction  between  the 
value  of  b and  the  physical  threshold  level.  However 
the  local  averaging  performed  by  the  LAS  method 
results  in  a slight  decrease  of  the  variance  of  the  dis- 
cretized field,  as  dictated  by  local  averaging  theory. 
In  that  each  field  is  represented  as  a discrete  lattice 
of  128  X 128  “cells,”  the  variance  of  each  cell  varies 
from  0.971  at  the  smallest  scale  of  fluctuation  con- 
sidered to  0.999  at  the  largest.  The  distinction 


where  <P  is  the  standard  normal  distribution  func- 
tion. The  estimate  of  E[Db],  denoted  mo*,  derived 
using  the  simulation  results  is  shown  in  Fig.  2 and  is 
in  complete  agreement  with  Eq.  (6).  Figure  2 (b) 
shows  the  estimated  standard  deviation  of  Db  de- 
noted sob-  Note  that  while  E[Db]  is  independent  of 
6,  its  variance  is  not.  In  keeping  with  the  practice  of 
normalizing  all  results,  the  scale  of  fluctuation  has 
been  normalized  with  respect  to  L = in  these 
plots.  Note  also  that  the  horizontal  threshold  axis 
values  decrease  to  the  right  — this  is  because  the 
threshold  axis  is  associated  with  time  in  the  next 
section  and  time  increases  to  the  right  as  usual. 

Figure  3 shows  the  estimated  mean  and  standard 
deviation  of  the  number  of  isolated  excursion  re- 
gions, Nb,  denoted  and  sn,,  respectively,  as  a 
function  of  scale  and  threshold.  Although  not  read- 
ily apparent  from  the  plot,  the  limiting  value  of  niNh 
at  the  right  edge  of  the  plot  (6^  — <»  ) is  1.  In  other 
words,  a single  excursion  exists  over  the  entire  do- 
main for  very  low  thresholds.  At  the  other  extreme 
(fo^oo),  the  limiting  value  is  zero  as  expected. 
Although  the  roughness  of  the  estimate  SNb,  as  seen 
in  Fig.  3(b),  is  as  yet  unexplained,  it  appears  un- 
likely that  it  arises  from  statistical  uncertainty  in  the 
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Fig.  2.  Estimated  statisties  of  total  excursion  area  per  unit  area,  D*,:  a)  mean,  b)  standard  deviation. 


Fig.  3.  Estimated  statistics  of  number  of  isolated  excursions,  A*;  a)  mean,  b)  standard  deviation. 


estimation  procedure.  Based  on  a sample  size  of 
2000  realizations,  a 90  percent  confidence  bound  on 
the  standard  deviation  of  Nt,  is  about  ±0.025^^* 
which  is  smaller  than  the  size  of  most  of  the 
“bumps”  seen  in  Fig.  3(b). 

Within  a given  realization,  the  average  area  of  an 
isolated  excursion  per  unit  domain  area,  De,  can  be 
obtained  using  the  number  of  excursions. 


4(b)  is  derived  using  the  assumption  that  the  sizes 
of  isolated  excursions  are  independent.  While  this  is 
true  from  realization  to  realization,  it  is  clearly  not 
true  within  a single  realization.  Thus  Fig.  4(b)  can 
only  be  considered  to  be  a rough  indication  of  the 
true  variability  of  the  area  of  isolated  excursions. 

3.  Calibrating  the  Model 


De  = 


Mb' 


Since  Db  is  the  sum  of  the  Nb  isolated  excursion 
areas,  the  expected  value  of  Dt  is  just 


£[De] 


E\Db} 


Figure  4(a)  illustrates  this  result  using  the  esti- 
mated mean  value  of  Nb  shown  in  Fig.  3(a).  The 
estimated  standard  deviation  of  De  shown  in  Fig. 


There  are  essentially  only  two  parameters  in  the 
random  field  model  considered  in  this  study.  These 
are  the  scale  of  fluctuation,  6,  and  the  threshold 
level,  h.  Although  there  is  little  published  experi- 
mental evidence  to  allow  a clear  statement  of  what 
these  parameters  should  be  for  a given  structure, 
some  preliminary  estimates  are  possible.  First  the 
relationship  between  the  threshold  level  and  the 
age  of  the  structure  can  be  obtained  by  fitting 
delamination  versus  age  data  collected  by  Trow  [7] 
and  presented  in  Ref.  [3]  to  the  curve  of  Fig.  2(a). 
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A good  fit  was  obtained  using  least  squares  re- 
gression (r^  = 0.93)  on  the  linear  relationship 
b=  3.17-0.225t.  To  obtain  this  relationship,  the 
total  fraction  of  delaminated  area  observed  by 
Trow,  Db,  was  plotted  against  b ~ 0~\1-Db)  and 
a best  fit  line  obtained.  Figure  5(a)  shows  the  re- 
gression results  and  Fig.  5(b)  shows  where  the  ob- 
served delamination  results  would  appear  on  Fig.  2. 
To  expand  the  scale,  only  positive  thresholds  are 
shown  in  Fig.  5(b).  The  choice  of  function  relating 
time  and  threshold  is  largely  arbitrary  as  long  as 
b{t)  is  a decreasing  function  over  all  times  of  inter- 
est (assuming  delaminated  areas  cannot  “heal”). 
While  a quadratic  gives  a slightly  better  fit  to  the 
raw  data,  it  violates  this  principle  and  so  cannot  be 
used  to  extrapolate. 


Some  additional  unpublished  data  was  made 
available  to  the  author  by  Public  Works  Canada 
based  on  a survey  of  a single  parking  garage  struc- 
ture in  Ontario,  Canada.  Figure  6 illustrates  this 
data.  The  ± 1 standard  deviation  curves  are  ob- 
tained using  Fig.  2(b)  with  OIL  =0.1.  Clearly,  for 
such  a choice  in  6,  the  variability  in  Db  predicted  by 
the  random  field  excursion  model  underestimates 
the  variability  in  the  observations.  The  additional 
variance  arises  because  the  time-threshold  relation- 
ship b{t)  is  itself  a function  of  random  coefficients, 
as  implied  by  the  regression  analysis.  Alternatively 
and  equivalently,  the  additional  variability  could  be 
ascribed  to  the  fact  that  the  delamination  field  is 
not  homogeneous  on  such  a scale.  When  considered 
at  the  scale  of  a typical  bay,  some  bays  show  much 


Fig.  4.  Estimated  statistics  of  isolated  excursion  areas  per  unit  area,  D^,:  a)  mean,  b)  standard  deviation. 


Fig.  5.  Regression  results  using  Trow’s  delamination  rate  data:  a)  regression,  b)  observed  versus  predicted  D*. 
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higher  delamination  rates  than  others,  perhaps  cor- 
responding to  smaller  mean  cover  depths  and/or 
higher  Cl~  inputs.  Note  also  that  larger  values  of  0/ 
L lead  to  higher  standard  deviations  in  Db  so  that 
the  variability  seen  in  Fig.  6 could  also  be  explained 
by  larger  values  of  6.  However  the  choice  of  6 
should  not  be  based  on  the  variability  in  Db,  rather 
this  observation  indicates  that  perhaps  the  Gauss- 
Markov  covariance  model  is  not  appropriate. 

In  this  study,  b(r)  will  be  treated  as  a determinis- 
tic function  which  corresponds  to  the  choice  of  a 
homogeneous  random  field.  It  is  beyond  the  scope 
of  this  initial  investigation  to  consider  estimating 
the  parameters  of  a nonhomogeneous  (or  self- 
similar) field  to  represent  the  delamination  process 


generally,  although  this  appears  indicated.  Since 
plan  views  of  the  delaminations  are  available  from 
Public  Works  for  individual  bays  of  approximately 
5 m X 5 m in  extent  (L  = 5),  attention  will  be  re- 
stricted to  data  on  such  a scale  over  which  the  field 
can  be  considered  homogeneous.  Choosing  one 
such  area  from  the  Public  Works  data,  the  mea- 
sured delamination  fraction  as  a function  of  time  is 
shown  in  Fig.  7.  In  this  case  the  line  of  best  fit  was 
found  to  be 

b=4.13-0.172r  (7) 

with  r^  = 0.86.  Notice  that  for  such  a case,  the  vari- 
ability in  the  observations  is  substantially  reduced 


Fig.  6.  Delamination  rate  data  provided  by  Public  Works  Canada.  The  fitted  line  is  b =3.61-0.142/  with  r^  = 0.6. 


Fig.  7.  Dclimination  rate  provided  by  Public  Works  Canada  over  a single  slab  region. 
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and  that  2/3  of  the  observations  lie  on  or  within  the 
±1  standard  deviation  curves  for  0/L=O.l  (see 
also  Fig.  5).  This  result  appears  encouraging  al- 
though it  is  recognized  that  it  could  be  due  in  part 
to  the  reduced  number  of  samples,  even  though 
similar  results  were  found  for  most  other  bays. 

Before  considering  the  estimation  of  9,  it  is  worth 
pointing  out  a further  difference  between  the  ran- 
dom excursion  model  and  the  commonly  accepted 
delamination  model.  Salt  induced  reinforcement 
corrosion  is  generally  believed  to  involve  two  stages 
(see  Ref.  [8]):  1)  an  initiation  phase  during  which 
the  alkalinity  of  the  concrete  surrounding  the  rein- 
forcement (which  renders  the  steel  passive)  is  re- 
duced by  the  migrating  Cl“  ions,  and  2)  an  active 
phase  during  which  corrosion  takes  place.  Attwood 
et  al.  [3]  estimates  the  initial  phase  to  last  4.7  years. 
During  this  phase  no  corrosion  is  assumed  to  take 
place.  Note  that  such  a model  can  only  be  applied 
to  points  in  the  slab  where  the  concrete  is  in  contact 
with  the  reinforcement. 

In  contrast  with  this  two  stage  model,  the  random 
field  excursion  model  admits  some  probability  of 
delamination  even  at  time  t =0.  Using  Eq.  (7),  one 
obtains  b =4.13  at  time  t =0,  so  that  the  expected 
total  delamination  area  per  unit  area  is  1.7  x 10“^  or 
about  17  mmVm^  of  slab.  At  these  levels  the  precise 
definition  of  delamination  comes  into  question.  If  it 
is  strictly  interpreted  as  a loss  of  bond  between  the 
concrete  and  reinforcement  then  this  result  is  not 
unreasonable  given  the  presence  of  cracks,  voids 
and  the  initial  state  of  the  reinforcement.  At  time 
t~5  years,  the  expected  total  delamination  area  is 
still  only  about  6 cmVm^  of  slab,  a level  which  is 
probably  still  largely  undetectable  at  the  surface  of 
the  slab  and  presumably  would  correspond  to  corro- 
sion in  the  immediate  neighborhood  of  surface 
cracks.  Although  it  is  believed  that  cracks  do  not 
contribute  significantly  to  the  areal  delamination 
process  [8],  it  is  not  unreasonable  to  expect  that 
they  can  be  initiators  of  the  corrosion  process  at  dis- 
crete points  in  the  slab.  If  this  is  the  case,  then  the 
many  year  delay  before  the  onset  of  observable  lev- 
els of  delamination  implies  that  the  corrosion 
growth  should  be  quite  slow  at  first.  The  results 
predicted  by  Fig.  2(a)  are  in  basic  agreement  with 
this  in  that  F[Dfc]  grows  very  slowly  for /)  decreasing 
to  about  2 [corresponding  to  t < 12  years  using 
Eq.  (7)]. 

Turning  now  to  the  estimation  of  the  scale  of 
fluctuation,  6,  it  becomes  apparent  that  this  task  is 
complicated  by  the  type  and  quality  of  data  avail- 
able. Ideally,  one  would  take  measurements  of  the 
corrosion  induced  stress  field  over  a number  of 


structures,  estimate  a spatial  covariance  structure 
and  from  this  obtain  6.  Even  if  this  approach  were 
possible,  the  nonhomogenieties  mentioned  above 
would  make  it  difficult.  However,  in  general  the 
stress  field  is  unmeasurable  and  what  little  data  is 
available  generally  consists  of  surveys  giving  the 
spatial  extent  of  delaminated  regions.  Figure  8(a)  il- 
lustrates such  a survey  while  Fig.  8(b)  is  a realiza- 
tion of  the  random  field  excursions  using  0 = 0.5  m. 
Once  the  time-threshold  relationship  has  been 
established,  and  for  the  purposes  of  this  argument 
Eq.  (7)  will  be  used,  a possible  technique  of  estimat- 
ing 0 would  be  to  count  the  average  number  of 
excursions  and  enter  Fig.  3(a)  at  the  appropriate 
threshold  to  estimate  OIL.  Purely  on  the  basis  of 
Fig.  8(a)  this  yields  estimates  of  0 = 2 to  4 m (0/ 
L =0.4  to  0.8).  However  realizations  at  this  scale 
yields  excursions  which  are  generally  far  too  well 
connected  as  shown  in  Fig.  9.  Realizations  at  such 
large  scales  appear  like  large  land  masses  with  many 
small  islands  close  to  shore.  Figure  8(a)  has 
“islands”  that  are  more  uniform  in  size  and  distri- 
bution implying  a smaller  scale  of  fluctuation.  In 
Fig.  8(b),  produced  using  0 = 0.5,  the  larger  islands 
are  of  similar  size  to  the  delaminations  seen  in  Fig. 
8(a).  This  along  with  arguments  to  follow  supports 
the  choice  of  a smaller  scale  of  fluctuation. 

While  one  cannot  expect  the  plots  in  Fig.  8 to  be 
identical  since  they  are  both  independent  realiza- 
tions of  a random  process,  any  more  than  one  could 
expect  the  pattern  of  delaminations  in  another 
building  to  be  identical,  a number  of  points  can  be 
made  about  the  two  plots; 

1)  there  is  no  apparent  spatial  orientation  of  the 
delamination  regions  in  the  observations  of  Fig. 
8(a),  indicating  that  the  assumption  of  isotropy  is 
acceptable,  at  least  for  this  case. 

2)  about  half  of  the  excursions  in  the  random  field 
model  (Fig.  8b)  are  of  very  small  extent.  This  frac- 
tion increases  at  larger  scales.  On  the  other  hand, 
in  Fig.  8(a)  there  are  only  very  few  “small’  areas 
appearing  in  the  later  survey.  It  seems  reasonable  to 
suspect  that  additional  small  delamination  regions 
are  in  fact  occurring  in  the  real  slab  but  that  the 
chain-drag  surveying  technique  is  unable  to  resolve 
them.  Operator  bias  will  almost  certainly  also  be 
present  due  to  the  prior  knowledge  of  existing  de- 
lamination areas. 

3)  the  random  field  model  is  much  “rougher” 
than  the  observed  delamination  plots.  Again  this  is 
likely  a problem  with  the  ability  of  the  chain-drag 
survey  to  resolve  detail. 
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Fig.  8.  a)  Observed  delamination  regions  at  ages  t = 13  years  (dark  grey)  and  t = 22  years  (light  grey)  on  a 
5mx5m  pwrtion  of  a parking  garage  slab,  b)  Excursions  of  a random  field  above  = 0.34(/ = 13)  and 
b = \9{t  = 22)  using  0 = 0.5  and  Eq.  (7). 


Fig.  9.  A realization  of  excursions  at  t = 13  years  (dark  grey) 
and  at  t = 22  years  (light  grey)  using  0 = 4 m on  a 5 m x 5 m field. 


The  last  two  points  illustrates  the  difficulty  in  esti- 
mating 6 on  the  basis  of  delamination  surveys.  The 
chain-drag  method  depends  on  setting  up  reverber- 
ations in  the  delaminated  concrete.  For  delamina- 
tion details  below  a certain  size,  the  frequency  shifts 
are  undetectable  to  even  the  most  sensitive  human 
ears,  rendering  them  unnoticeable.  Thus  many  of 
the  delamination  details  on  which  an  estimation  of 
0 depend  are  unavailable  using  current  surveying 


techniques.  To  some  extent  the  deficiencies  in  the 
chain-drag  procedure  could  be  accommodated  sim- 
ply by  introducing  more  local  averaging  in  the  ran- 
dom field  model  — effectively  smoothing  the  field. 
While  such  a “correction”  would  not  likely  result  in 
an  improved  delamination  model,  it  may  allow 
improved  estimates  of  0.  For  the  purposes  of  this 
investigation,  the  realization  based  estimate  of 
0 = 0.5  m is  used. 

4.  Discussion 

Conceptually  the  model  proposed  herein  is  quite 
attractive  in  that  the  delamination  process  is  indeed 
a threshold  excursion  process  in  two  dimensions. 
Once  the  details  of  the  model  have  been  established 
(type  of  distribution,  time-threshold  relationship 
and  scale  of  fluctuation)  and  some  of  the  properties 
of  threshold  excursions  in  two  dimensions  have 
been  determined  analytically  or  via  simulation,  the 
model  can  be  used  in  a reliability  context.  For  ex- 
ample, using  Eq.  (7),  at  time  t = 12  years  {b  =2.07) 
Fig.  2 along  with  Eq.  (6)  indicates  that  =0.0194 
and  sot  — 0.0080,  using  0 = 0.5  m.  If  Dt  is  assumed  to 
have  a Beta  distribution,  which  is  properly  bounded 
between  0 and  1,  then  the  probability  that  Db  is  less 
than  3 % at  4 years  is  given  by 

<0.03]  = 0.90  . (8) 

In  turn,  if  a target  reliability  of  0.9  were  chosen 
against  delamination  in  excess  of  3 %,  then 
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inspection  of  the  garage  would  be  recommended  at 
t = 12  years  for  this  structure. 

Because  the  model  includes  the  spatial  aspects  of 
delamination,  it  can  be  used  to  evaluate  testing  and 
mapping  procedures.  For  example,  at  time  r = 12 
years  one  could  expect  about  42  isolated  excursions 
on  a 5 mx5  m slab  (0  = 0.5  m)  with  a standard  de- 
viation of  about  11.  The  area  of  each  isolated  excur- 
sion would  average  about  0.011  m^  with  a standard 
deviation  of  about  0.022  m^.  Surveys  which  yield  re- 
sults considerably  different  than  these  may  require 
verification  using  other  techniques.  In  addition,  the 
spatial  description  of  delamination  could  be  used  in 
a structural  reliability  study.  For  this  the  simulation 
approach  could  be  employed  to  yield  a measure  of 
the  degree  of  clustering  of  the  delamination  regions 
(see  Ref.  [6]). 

One  recognizes  that  the  excursion  model  is  at- 
tempting to  predict  the  cracked  state  of  a concrete 
slab.  In  that  internal  cracks  are  exceedingly  difficult 
to  map,  even  in  controlled  laboratory  conditions, 
the  model  is  to  some  extent  intuitive  and  will  likely 
remain  so  until  improved  surveying  techniques  are 
developed.  Nevertheless  the  model  demonstrates 
some  promising  features  and  can  be  used  as  a pow- 
erful reliability  tool  when  its  parameters  are  clearly 
defined  in  terms  of  additional  data  and  simulation 
studies.  In  particular,  the  data  shown  in  Fig.  6 indi- 
cates that  perhaps  alternative  correlation  functions 
[see  Eq.  (2)]  should  be  studied  — multiple-scale  or 
self-similar  type  random  fields  are  suggested,  show- 
ing the  smalt  scale  behaviour  over  small  regions 
while  reflecting  also  the  large  scale,  slower  varia- 
tions over  larger  domains.  Also  the  fact  that  the 
model  allows  delamination  to  occur  at  time  r = 0 (or 
before)  implies  that  some  thought  should  be  given 
to  the  assumption  of  a Gaussian  random  field  and/ 
or  the  time-threshold  relationship.  In  the  interim, 
however,  the  choice  of  a Gaussian  field  and  linear 
time-threshold  relationship  leads  to  results  which 
appear  reasonable  for  any  time  t > 0. 
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